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Abstract. We report results of ^^As nuclear magnetic resonance (NMR) and nuclear quadrupole resonance 
(NQR) experiments as well as ^Li NMR on different samples of self flux grown LiFeAs and 5% Co doped 
LiFeAs single crystals, and a polycrystalline LiFeAs sample. We were able to distinguish the samples by 
their slightly different quadrupole frequencies, uq, which is a direct measure of the electric field gradient 
(EFG) at the As site. Interestingly, samples with a large quadrupole frequency appear to show a different 
Knight shift and spin lattice relaxation in the superconducting state from those with a lower uq, yet all 
the samples are clearly superconducting. For sample SI which has the largest uq, we find constant Knight 
shift /C across Tc for a certain direction of the magnetic field and a peculiar upturn of the NQR spin lattice 
relaxation rate (TiT)"^ below Tc- In contrast, samples with a lower i/q exhibit the expected behavior for 
a singlet superconductor: a drop of /C and {TiT)~^ for both NMR and NQR below Tc- Our results show 
that already tiny changes in stoichiometry uncovered by slightly different NQR frequencies lead to very 
different behavior of the NMR parameters in the superconducting state of LiFeAs. Different possibilities 
will be discussed which may explain the contrasting behavior. 

PACS. XX.XX.XX No PACS code given 



1 Introduction 

Among the recently discovered superconducting iron pnic- 
tides [HH], LiFeAs is a rare member exhibiting supercon- 
ductivity with Tc ^ 18 K in a stoichiometric form without 
dopants or pressure. Yet, the Li content in this compound 
is difficult to control and seems to have a large influence 
on the physical properties. For example, Li deficiencies of 
only ^1% greatly suppress superconductivity [3]. Proba- 
bly related to the sensitivity of the Li concentration, the 
superconducting ground state of LiFeAs and the pairing 
mechanism are still under debate. 

Small angle neutron scattering (SANS) and angle re- 
solved photoemission spectroscopy (ARPES) [4l[5l|6] sug- 
gest that LiFeAs could be a weakly electron-phonon cou- 
pled conventional-type superconductor. However, a recent 
Raman scattering study did not find evidence for substan- 
tial electron-phonon-coupling and no superconductivity- 
induced phonon anomalies [7]. Furthermore, while the spin 
density wave (SDW) state is absent in LiFeAs, there is ev- 
idence that weak local moments [5] and magnetic fluctua- 
tions are still present in the normal state, putting LiFeAs 
close to a magnetic instability [^fTU] . Theoretical analyses 
of the electronic band-structure find a superconducting 
order parameter of s± type driven by collinear antiferro- 
magnetic fluctuations [TT]. The fully gapped s± supercon- 
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ductivity is also suggested by heat transport [T^, mag- 
netic penetration depth measurements jl3j . and inelastic 
neutron scattering (INS) study [13]. On the other hand, 
Brydon et al. [TS] proposed spin-triplet p-wave pairing in 
this system. Being in agreement with the theoretical ar- 
gument, vortex properties jl6j and Hc2 measurements |17j 
for H \\ ab show a very similar behavior as in the suppos- 
edly triplet superconductor Sr2Ru04, and quasi particle 
interference (QPI) patterns measured by scanning tunnel- 
ing microscopy (STM) [l^ supports an elementary p-wave 
symmetry, rather than singlet pairing symmetries (s^- or 
d-wave) . 

The so far reported NMR measurements of LiFeAs 
powder samples indicate the existence of magnetic correla- 
tions from the analysis of the Korringa relation [TOll^ . In 
the superconducting state, the Knight shift shows a sharp 
drop at Tc which is suggestive of spin-singlet superconduc- 
tivity. In the only NMR study on single crystals Ma et al. 
argue that the absence of the spin density wave ordering in 
LiFeAs is due to ofF-stoichiometry and/or lattice defects, 
and find two different Li sites in their superconducting 
samples pi] . 

Here, we report detailed NMR and NQR results on 
three single crystal LiFeAs samples, as well as on a 5% Co 
doped single crystal, and a polycrystalline LiFeAs sam- 
ple. We show that the NQR frequency is highly sensi- 
tive to the Li content, and that samples with tiny differ- 
ences in stoichiometry can be distinguished by their differ- 
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ent quadrupolc frequencies. A doping dependence of the 
quadrupole frequency has already been found in other iron 
pnictides |22| , and similar changes of vq with Co doping in 
LiFeAs indicate that also Li deficiencies change the dop- 
ing level of the samples. Interestingly, samples with a large 
quadrupole frequency appear to show a different Knight 
shift and spin lattice relaxation in the superconducting 
state from those with a lower uq. Sample SI which has 
the largest vq, exhibits a constant Knight shift K. across 
Tc for a certain direction of the magnetic field and a pe- 
culiar upturn of the spin lattice relaxation rate (TiT)"^ 
below Tc- In contrast, samples with a lower vq exhibit a 
drop of /C and {TiT)~^ below T^- Our results show that al- 
ready tiny changes in stoichiometry uncovered by slightly 
different NQR frequencies lead to very different behavior 
of the NMR parameters in the superconducting state of 
LiFeAs. The tiny differences in stoichiometry and their 
large impact on the superconducting properties may also 
account for the contradicting results reported so far. 

2 Sample preparation and Experimental 
details 

The LiFeAs single crystals were grown by a self-flux method 
using a molar ratio of Li:Fe:As = 3:2:3 similar to ref. ^5] . 
The stoichiometry of the samples has been checked by in- 
ductively coupled plasma mass spectroscopy (ICPMS). 20 
mg of the LiFeAs single crystal were dissolved in a leak free 
glass ampoule in nitric acid. The molar ratio Li:Fe:As is 
found to be 0.99:1.00:1.00, consistent with a stoichiomet- 
ric LiFeAs composition [53]. ARPES measurements per- 
formed on a LiFeAs single crystal from the same batch 
were found to be in agreement with an exact stoichiome- 

try 0. 

The susceptibility of our samples exhibits a sharp su- 
perconducting transition at Tc ^ 18 K [^H]. A clear 
anomaly at Tc in the specific heat measurement [23], and 
the large residual resistance ratio (RRR) [231125] , with the 
lowest residual resistivity observed so far in iron pnic- 
tides [^17711^ indicate that the samples are clean su- 
perconductors. Furthermore, STM measurements yield a 
defect concentration < 1 % [TH] . Another method to check 
the sample quality is the measure of the linewidth of a 
NQR spectrum or, equivalently, of the satellite transitions 
{I^ = -3/2 o -1/2 and h = +1/2 ^ +3/2) of a NMR 
spectrum, because those measure directly the distribution 
of the electric field gradient (EFG) which is caused by 
disorder or defects in the sample, while the broadening 
of the central transition (/^ = —1/2 o +1/2) is domi- 
nantly magnetic in origin. As we will show below, all sam- 
ples exhibit very narrow linewidth of both As NQR lines 
and ^Li NMR satellites which prove the absence of large 
amounts of defects in our samples and a stoichiometry of 
close to 1:1:1. 

^5 As (nuclear spin 7 = 3/2) NMR/NQR, and ^Li (/ ^ 
3/2) NMR measurements were carried out in five differ- 
ent LiFeAs samples: four single crystalline samples: SI, S2, 
S3, and Co 5%-doped LiFeAs, and a polycrystalline sam- 
ple, where the single crystals SI, S2, and S3 are from the 



same batch. Due to the extreme sensitivity of the samples 
to air and moisture, all the samples were carefully sealed 
into quartz tubes filled with Ar gas. It could happen that 
thermal cycling damages the sealing of the quartz tubes 
that contain the sample, so that, if the sample probe has to 
be taken out of the cryostat, the sample is easily degraded 
by contact with air. For this reason we could not obtain 
full data sets for all of the samples since the NMR and 
NQR signals become negligibly weak in a sample which 
had contact with air. 

Since the local symmetry at the ^^As is axial (tetrag- 
onal), the nuclear quadrupole frequency vq can be deter- 
mined directly from the resonance frequency of the NQR 
spectrum for the ''^ As and by the splitting of the satellites 
for the ^Li. The Knight shift was obtained by measuring 
the central transition. The nuclear spin-lattice relaxation 
rates T^^ were measured by saturation and inversion re- 
covery methods. For exact orientation of the crystals with 
respect to the static magnetic field a single axis goniome- 
ter has been used for most of the measurements. 



3 Experimental Results 

3.1 ^5 As NQR 

For the ^^As in a non-cubic environment, the four-fold de- 
generacy of the nuclear spin of / — 3/2 is partially lifted by 
the interaction between the nuclear quadrupole moment 
Q and the surrounding EFG, eg, which allows a NQR res- 
onance at a frequency given by vq = e^qQjh where h is 
Planck's constant and e the electron charge. We find a very 
narrow NQR line near 21.5 MHz at room temperature, as 
shown in the inset of Fig. 1. The ^^As NQR spectra have 
a width ranging from 60 kHz (S2) to 80 kHz (SI) at room 
temperature, which is a factor of 2-3 narrower than 170 
kHz reported earlier in powder samples of LiFeAs |20j . 
and even significantly smaller than the values reported 
for other undoped (non-superconducting) iron pnictides: 
100 kHz from NMR satellites in BaFe2As2 [H], 308 kHz 
and 385 kHz in La- and SmOFeAs, respectively [35], and 
480 kHz in CaFe2As2 [SO]. This is even more surprising 
regarding the high quadrupole frequency of LiFeAs with 
respect to the other iron pnictides, and indicates the high 
homogeneity of all of our single crystals. 

The temperature dependences of the nuclear quadrupole 
frequency vq of the ^^As are shown in Fig. 1. The re- 
sults obtained in the polycrystalline sample and the 5% 
Co doped single crystal sample are also compared. We find 
that VQ^ or the EFG eg, varies with the samples investi- 
gated. Despite the different vq at a given temperature, 
the temperature dependences of vq of all the samples are 
very similar among one another, i.e., vq decreases with 
decreasing T, being saturated at low temperatures. Such 
a strong temperature dependence of vq is commonly ob- 
served in iron pnictides [291I31| , and is attributed to mostly 
electronic effects due to complicated multi-orbitals of the 
iron ion (3(i^), since the lattice (thermal vibration) con- 
tribution, usually, leads to an opposite (and very weak) 
temperature dependence. However, we can distinguish the 
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Fig. 1. The quadrupole frequency, vq, decreases with decreas- 
ing temperature and almost saturates at low temperatures. 
While the temperature dependences of of all the samples 
are similar, the values of the sample SI are noticeably larger 
than those of the samples S2 and S3, as clearly shown in the 
inset. For comparison, vq{T) for poly crystalline and Co-doped 
samples are shown, too. Note that the linewidth is also strongly 
sample-dependent, yielding 80, 60, and 40 kHz for SI, S2, and 
S3, respectively at room temperature. 



samples clearly by their different values at a given tem- 
perature, which resembles a doping dependence similar 
to other pnictides even though samples S1-S3 are 

from the same batch, and should have the same compo- 
sition according to ICPMP and ARPES [HIS]. Interest- 
ingly, we observe that vq (T) of the polycrystalline sample 
is very similar to that of SI, while vq{T) of the Co-doped 
sample is even below those of S2 and S3. It appears that 
samples with a larger vq reveal a peculiar temperature de- 
pendence of the relaxation rates and/or the Knight shift 
in the superconducting state, as shown below, whereas 
samples with a lower vq show more normal behaviors, as 
expected in typical spin-singlet superconductors. We note 
that yqiT) of the powder samples obtained by Li et al. 
[20] appears to be close to our Co-5% doped single crystal. 

Fig. 2 shows the ^^As NQR measurement of the spin 
lattice relaxation rate divided by temperature, {TiT)~^. 
For all samples, (TiT)~^ slightly decreases with decreas- 
ing T and approaches a constant value below ~ 150 K 
down to Tc- Below Tc, however, the temperature depen- 
dence of (TiT)~^ varies among the samples. Due to the 
opening of the superconducting gap, the electron density 
of states at the Fermi level to which (TiT)~^ is propor- 
tional, is reduced rapidly, and thus {TiT)~^ is expected 
to decrease accordingly. While S3 exhibits such a rapid 
drop at Tc, SI does not. Nevertheless, the fact that the 
upturn occurs at Tc suggests that the unusual behavior 
is associated with superconductivity, as will be discussed 
below in comparison to the NMR (TiT)"^. For the Co- 
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Fig. 2. ^^As NQR spin-lattice relaxation rates divided by T, 
(riT)~\ as a function of T in zero field. While (TiT)"^ of all 
four samples are essentially the same above Tc, those change 
drastically below Tc. For both the single crystal SI and the 
polycrystalline sample, {T\T)~^ rises below Tc. 5% Co-doping 
suppresses the enhancement below Tc, but still without a de- 
crease. On the contrary, data from another single crystal S3 
reveals a rapid drop, as expected in the superconducting state. 



doped sample, the enhancement of (TiT)^^ is suppressed, 
while maintaining the constant value from above T^. The 
polycrystalline sample also shows a strong enhancement of 
(TiT)"^ below Tc, similar to that of SI. In all cases, su- 
perconductivity has been confirmed at the same time by 
a strong change in the resonance frequency of the NMR 
circuit, which is proportional to the ac susceptibility of 
the sample. We note that both SI and the polycrystalline 
sample which reveal a similar behavior of (TiT)~^ possess 
larger vq values than the other samples, as shown in Fig. 
1. Tentatively, this trend may suggest that the strength 
of the EFG has a relationship with the underlying physics 
which may cause the different behaviors of (TiT)~^ in the 
superconducting state. 



3.2 ^ Li NMR spectra 

Fig. 3 shows ^Li NMR spectra including the two satellites 
obtained for the crystal S3. At 200 K, the spectra for both 
directions along c and ah reveal very sharp lines which al- 
low us to determine vq = 33 kHz accurately. This value is 
also consistent with the value of vq = 34 kHz estimated 
by Jeglic et al. [19j by an echo decay measurement. The 
almost identical linewidth of cental and satellite lines indi- 
cates that the broadening mechanism is mainly magnetic, 
and the distribution of the EFG is negligibly small. At low 
temperatures, the ^Li spectra broaden significantly, unlike 
the ^^As NMR-NQR spectrum whose linewidth increases 
only for SI at low temperatures. This indicates that the 
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Fig. 3. ^Li NMR spectra at 200 K and 20 K at H = 4.9994 
T, witii tiie ^Li quadrupole frequency, uq — 33 kHz. Tiie sep- 
aration between the satellite lines for H \\ c corresponds to 
2i'Q, as expected in the tetragonal symmetry (the asymmetry 
parameter rj = 0). The linewidth, which is the same for central 
and satellite lines, for H \\ ah {H \\ c) increases from 9 kHz (11 
kHz) at 200 K to 20 kHz (40 kHz) at 20 K. 



^Li nuclei experience quasi-static spin fluctuations at low 
temperatures, particularly for H along the c axis. 

We emphasize that our ^Li NMR spectra assure that 
there is only one single Li site in our crystals. This is 
in stark contrast with the crystals used by Ma et al. [H] 
where two ''Li resonances are observed, indicating the pres- 
ence of inequivalent Li sites in their crystals. Therefore our 
results are not consistent with the claim that the absence 
of magnetism is due to the ofF-stoichiometry or defects. 
Rather, we argue in this paper that the stoichiometry is a 
critical parameter governing the nature of magnetism and 
superconductivity in LiFcAs. 



3.3 ^^^As NMR Knight shifts 

Fig. 4 (a) shows the T-cvolution of the ""^As NMR central 
line of the single-crystal SI measured with H parallel to 
the a6-plane. At high temperatures, the spectrum shifts 
down with decreasing T, and approaches a constant fre- 
quency below ^ 50 K. Unexpectedly, we observed that the 
spectrum (drawn with thick lines below Tc) maintains the 
resonance frequency in the superconducting state down 
to 4.2 K. This behavior drastically changes when H is ro- 
tated out of the afc-plane by ~ 2° - now the As NMR line 
shows visible shifts below T^, see Fig. 4 (b). Nevertheless, 
this indicates that the sample is superconducting. 

Due to the large NQR frequency vq, there should be 
a substantial shift of the central transition for H \\ ab 
that is strongly angle dependent, due to the second order 
quadrupole effect. In an uniaxial symmetry, which we con- 
firmed by ^Li NMR spectra, the second order quadrupole 
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Fig. 4. ^'''As NMR spectra of the central transition (/^ — 
1/2 o -1/2) at an external field oi H ^ 7.0494 T as a func- 
tion of temperature for the single crystal SI. Spectra in the 
superconducting state are shown as thick lines, (a) Spectra for 
H II ab. Note that the resonance frequency of the spectra does 
not change upon going through Tc. (b) Same as in (a), but 
with tilting the sample by 2°. Although the temperature de- 
pendence of the shift above Tc is similar to the case of the 
exact alignment, the resonance frequency drops at Tc{H) ~ 13 
K, indicating a vanishing spin susceptibility. 



shift for a nucleus with spin / = 3/2 is given by |33j 



(1 - cos2 6i)(i - gcos^e*), (1) 



where 9 is the angle between the tetragonal c axis and H. 

Ai'{9) has been subtracted from the total shift of the 
NMR lines to extract the Knight shift shown in Fig. 5. In 
our case, the extreme sensitivity of Ai/ to 6 caused by the 
large vq is indeed a benefit since one can take advantage 
of it in order to align the sample with great accuracy. In 
particular, it is very useful for H \\ ab, since Av is simply 
the maximum when = 90°. Furthermore, we confirmed 
by exact diagonalization of the nuclear Hamiltonian the 
validity of the second order correction. 

The resulting Knight shift K. shown in Fig. 5(a) de- 
creases slowly upon lowering the temperature, flattening 
out at around 50 K, being in agreement with the results re- 
ported previously pOlfT^ . as already reflected in the raw 
data (see Fig. 4). The total Knight shift K. consists of 
a spin and an orbital contribution, /C = /Cspin + ^orb, 
where the latter involves the orbital motion of the con- 
duction electrons and is usually temperature independent. 
Here, M^spin is directly proportional to the spin suscepti- 
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bility, /Cspin = ^Xspin, where A is the hyperfine couphng 
strength. Therefore, /C should vanish at T ^ Tc if the 
spin state of Cooper pairing is singlet. Usually, the hyper- 
fine coupling constant is extracted from plots of /C versus 
the bulk susceptibility, x- We extracted the hyperfine cou- 
pling constants for temperatures T > 170 K, because be- 
low the temperature the Knight shift does not scale with 
the susceptibility x due to a magnetic impurity contribu- 
tion which affects only x- For H \\ ah we obtain Aab = 6.3 
T/hb, and for H \\ c Ac = 0.95 T/^b- These values ap- 
pear to be strongly anisotropic, compared to other iron 
pnictides, e.g., Aab = 3.87 T//x_b and Ac = 2.61 T/^s 
for isostructural NaFeAs [34] and Aab = 2.64 T/fis and 
Ac 1.88 T/iiB for BaFezAsz [Hj. 

A scaling of the Knight shifts measured for different 
orientations and for different nuclei would indicate that all 
nuclei probe the same component of the spin susceptibil- 
ity. Such a behavior has been found for optimally doped 
LaOo.gFo.iFeAs I35|. Instead, we find that for ^^As IC"'' 
does not give a linear relation with /C^Q Likewise, the 
Knight shift of the ^Li nucleus does not scale with the 
Knight shift measured at the ''^As for both directions. 
Since there are several bands crossing the Fermi level in 
the iron pnictides, one can expect each band to make a 
different contribution to the total spin susceptibility with 
different hyperfine couplings to different bands. Such is 
the case for ^^O NMR in Sr2Ru04, where the oxygen 
simultaneously couples to multiple bands with different 
temperature-dependent susceptibilities [351 ■ The strong an- 
gular dependence of /C is also observed for the ^^O Knight 
shift in Sr2Ru04, whose sign changes with angle at low 
temperatures as in our case. While the spin susceptibil- 
ity Xspin is always positive, the Knight shift can be nega- 
tive since /Cspin can be decomposed into two components: 
/Cspin = AsXs+AcpXnon-s whcrc Ag is the direct Fermi con- 
tact hyperfine coupling to s-electrons and ^cp arises from 
core polarization of inner s-shells due to non-s electrons 
{p or d) [37] ■ Here, As is always positive whereas ^cp is 
always negative. Thus, if Xnon-s is strongly angle and tem- 
perature dependent, K. will change accordingly, possibly 
reversing its sign. Although this may catch the essential 
features, the whole T- and angle-dependencies of our data 
obtained in LiFeAs are not correctly understood quantita- 
tively within this simple picture. This is likely due to the 
multi-band structure as mentioned above. In this case, the 
spin susceptibility from each band may have quite differ- 
ent response to temperature and field direction, depending 
on the overlap with p-orbitals of the As ion, resulting in 
the observed angular and temperature dependencies of the 
Knight shift. 

Now, we focus on the temperature dependence of /C 
at low temperatures, which is shown in Fig. 5(b). Upon 



^ We do not show IC^ here, because we measured our sam- 
ples only with a single axis goniometer. Therefore, while one 
can reach a perfect alignment for H \\ ab, there may be small 
deviations for H \\ c which lead to additional second order 
quadrupole correction. This affects only the absolute value of 
/C^, but not the temperature dependence which can be com- 
pared with /C''\ 



going through the superconducting Tc, the Knight shift 
K."'^ for the sample SI does not show any change, which 
was already manifested in the raw data [Fig. 4(a)]. 

Furthermore, in the superconducting state the behav- 
ior of K, changes when the field H is tilted by just 2° out 
of the ab plane: there is a clear drop of K. at Tc and it 
approaches to a finite value as T ^ (see Fig. 5). For an 
angle of 5.6°, JC is small and even upturns slightly below 
Tc, but approaches a similar value for 2°-off case, too, as 
T — > 0. The temperature dependence of /C for a tilting of 
2° is similar to the data measured in the polycrystalline 
sample [TOll^ . and further indicates that our crystals are 
superconducting. In polycrystalline samples, however, one 
may argue that whether K.""^ is constant for precisely in- 
plane field is very difficult to be confirmed due to the in- 
evitable ambiguity in assigning the singularity for H \\ ab 
in a powder pattern. 

However, the results measured in another single crystal 
S2 exhibit totally different behavior. IC for S2 is slightly 
larger than for SI above Tc, and /C°^ shows a clear drop at 
Tc- With tilting the crystal by 3°, a similar behavior of K, 
was observed, unlike the strong angle dependence obtained 
for SI. Such a strong variation of /C in two samples grown 
in the same batch is surprising, and we interpret this as 
an indication that LiFeAs is located in the vicinity of an 
instability which may cause the extreme sensitivity of the 
physical properties. 

3.4 ^5 As NMR and ISIQR linewidth 

Fig. 6 shows the temperature dependence of the full width 
at half maximum (FWHM) for the NMR and NQR spec- 
tra. For sample SI, FWHM increases with decreasing tem- 
perature, indicating a growing of magnetic fluctuations. In 
contrast, FWHM of sample S2 is nearly temperature inde- 
pendent. Note that for technical reasons we could obtain 
only limited high temperature NMR data for S2 (see also 
Section[2]). The T independent NQR FWHM supports the 
fact that the NMR FWHM is also T independent, espe- 
cially since its absolute value is very small even at low 
temperatures. 

In general, the sample quality is one of the factors de- 
termining the linewidth of NMR and NQR spectra. When 
interpreting the linewidth one has to distinguish intrinsic 
effects such as broadening from local magnetic moments 
from the effect of impurities. This can be done by com- 
paring the temperature dependence of the linewidth of 
NMR and NQR spectra where the broadening from de- 
fects (mainly quadrupole effects) and local magnetic mo- 
ments (magnetic effects) appear differently. In the para- 
magnetic phase, the local moments are fluctuating. If the 
inverse correlation time, I/tc of the fluctuations is much 
higher than the NMR linewidth, each nucleus sees only the 
time-averaged local field and the NMR line should be nar- 
row (motional narrowing). As the fiuctuations slow down, 
the nuclei will start to feel a distribution of local fields 
which broadens the NMR/NQR lines. Note that even a 
tiny moment size can broaden the NMR line consider- 
ably. For example, if we assume a moment size of only 
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Fig. 5. (a) Tlie Kniglit sliift (/C) as a function of temperature 
for sample SI for different orientations, and for S2 for H \\ ab. 
(b) JC at low temperatures. For SI and H \\ ab, K. is constant 
across Tc, whereas it decreases below Tc as soon as the field 
is tilted out of the ab plane. In contrast, for the crystal S2, K, 
drops below Tc even for H \\ ab, and the angle dependence of 
K, is not as strong as that of SI. 



lO^^^B, it will lead to a linewidth on the order of ^ 18 
kHz, from the relation Av ~ 7?int A^l, using the hy- 
perfine coupling A = 6.3 T/fJ.B as extracted in Section 
13.31 The temperature dependence of FWHM (Fig. 6) cor- 
roborates the magnetic broadening. Disorder or defects in 
the sample can also lead to an enhanced magnetic broad- 
ening, but the effect of disorder is much stronger on the 
quadrupole broadening, and would not lead to the ob- 
served temperature dependence. However, the influence 
of quadrupole effects can be directly measured by NQR 
or by the satellite transitions (/z = —3/2 o —1/2 and 
/z ^ +1/2 ^ +3/2) of the NMR spectrum, whereas the 
central transition (/j, = —1/2 -f^- -1-1/2) is only affected 
by second order quadrupole effects. Any small deviations 
from a homogeneous charge distribution or lattice anoma- 
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Fig. 6. FWHM as a function of temperature. For sample 
SI, FWHM measured by NMR and NQR increases with de- 
creasing temperature, and shows a strong angle dependence. 
In contrast, for the crystal S2, FWHM is smaller than for SI 
and the NQR FWHM is temperature independent. The NMR 
and NQR FWHM of S2 seem to be T independent as indicated 
by the dashed line. In the superconducting state, FWHM for 
SI decreases, whereas it increases for S2. 



lies such as deficiencies or defects will lead to a distri- 
bution of EFGs at the nucleus which contribute to the 
quadrupole linewidth. Therefore, while the narrow NQR 
lines of our LiFeAs single crystals is related with high ho- 
mogeneity, the increase of the linewidth when approaching 
Tc, regardless of the rotation angle, signals that magnetic 
correlations progressively gain in strength upon lowering 
the temperature in the sample SI, whereas S2 does not 
show such a behavior. 

We observe that the narrow line for H \\ ab rapidly 
broadens when H rotates out of the plane. A tilting an- 
gle of 2° already causes a noticeable broadening (Fig. 6) 
and for 5.6° the line broadens by a factor of three. For 
about 8° the width exceeds 200 kHz and K, is very small. 
Such an angle dependence could arise from a magnetic ef- 
fect, where anisotropic, probably momentum-dependent, 
spin fluctuations are present, which are closely related 
to the anisotropic Knight shift K,. This would be in line 
with the presence of incommensurate, nearly ferromag- 
netic, spin fluctuations that emerge in microscopic calcu- 
lations jT5| , the character of which can be strongly affected 
by a magnetic field. Or the angle dependence arises from 
a quadrupolar effect. In this case, a distribution of angles 
between the external magnetic field H and the direction of 
the principle axis of the EFG tensor Vzz would be needed 
to explain the angle dependent broadening, yet narrow 
NQR resonance line. 

In the superconducting state, the linewidth for sample 
SI unexpectedly decreases below Tc- This feature is also 
observed in the NQR spectra. Usually the NMR linewidth 
increases in the superconducting state due to vortex-related 
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broadening as is observed for sample S2. The origin of the 
decreasing linewidth just below Tc is not yet clear. Nev- 
ertheless, it indicates that the magnetic properties of the 
sample change at the superconducting transition temper- 
ature, namely that the spin fluctuations which lead to the 
broadening above Tc are strongly reduced in the super- 
conducting state. 



3.5 '^^As spin-lattice relaxation rates in external field 

By As NQR, we have shown that the ^'^As (TiT)~i is 
strongly sample dependent in the superconducting state. 
In order to check how an external field affects the low en- 
ergy spin dynamics in the superconducting state, we mea- 
sured {TiTy^ by ^5 As NMR in the two single crystals SI 
and S2. As shown in Fig. 7, for the crystal SI, we observe 
that (TiT)~^ even in the normal state is enhanced com- 
pared to the NQR results (solid horizontal line). Since it 
is expected that the NMR results for || c are equivalent 
to NQR ones for which the nuclei are also quantized along 
the c axis, the behavior is somewhat surprising, suggesting 
the presence of unusual, field dependent low energy spin 
dynamics in the system. For H \\ afo, (TiT)^^ is further 
enhanced, as denoted by the dotted line. Regardless of the 
direction of if, the strong enhancement of (TiT)"^ below 
Tc observed in the NQR measurement (i.e., in zero field) 
is substantially suppressed in a magnetic field, although a 
sharp drop is still absent for sample SI. A similar behav- 
ior has been found in Lao.syCao.iaFePO [35|, where '^^P 
{TiT)~^ increases just below Tc in low magnetic fields, 
while high magnetic fields of about 6 T suppress this in- 
crease (see also the discussion in Section!?]). 

For the crystal S2, we observe contrasting results, as in 
the Knight shift (Fig. 5). In the normal state (TiT)"! data 
for H \\ ah are considerably smaller than those of SI. Since 
NQR results were the same among all of the samples, this 
indicates that the spin dynamics even in the normal state 
become material-specific in field. In the superconducting 
state, a sharp drop at Tc is observed for S2, as in NQR 
case for S3. 

An interesting observation is that (TiT)^^ of SI re- 
veals a weak but clearly visible maximum just above Tc 
for both directions of H . Together with the enhanced nor- 
mal state (rir)^"'^ in an external field, one can argue that 
the normal state spin dynamics are strongly influenced 
by an external field, signaling a nearby critical instability 
whose nature is unclear yet. Furthermore, the presence of 
the local maximum of (TiT)~^ just above Tc in field as 
well as the mysterious upturn of (TiT)~^ below Tc in zero 
field may indicate unusual superconductivity in LiFeAs. 



4 Discussion 

Our single crystals SI, S2, and S3 of LiFeAs which were 
grown in the same environment, exhibit significantly dif- 
ferent static {uQ and /C) and dynamic [(TiT)~^] prop- 
erties, particularly, in the superconducting state. Never- 
theless, all the crystals arc found to be of good quality 
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Fig. 7. ■^"As NMR (TiT) 



T(K) 

as a function of T measured at 
7 T for the sample SI and at 8.5 T for the sample S2. Above 
Tc we observe a small but visible enhancement compared to 
NQR results (horizontal solid line), particularly for SI. Below 
Tc, only the sample S2 exhibits a drop of (TiT)"^. Sample S3 
exhibits similar temperature dependence as sample S2. 



and high homogeneity as evidenced by the narrow ^^As 
NQR line and the almost equal linewidth of central and 
satellite ^Li NMR lines. Ahhough the NQR FWHM of SI 
(-80 kHz) is much larger than that of S3 (-40 kHz), it 
is still a factor of more than 2 narrower than reported for 
powder LiFeAs samples [50] and significantly smaller than 
reported for other undoped iron pnictides (sec above). 
Therefore, our results suggest that pure LiFeAs is located 
near a critical point so that its physical properties are ex- 
tremely sensitive to very small perturbations such as tiny 
variations in stoichiomctry. A similar critical doping de- 
pendence has also been reported in literature [3], where 
one percent of Li deficiencies largely suppress supercon- 
ductivity. However, we emphasize that all of our crystals 
are superconducting with similar Tc- We confirmed the 
onset of superconductivity by a strong change in the reso- 
nance frequency of the NMR sample probe below Tc which 
is caused by the change of the surface resistance of the 
sample. Furthermore, a drop of the Knight shift for small 
angles off || ab and the anomalous change of both the 
linewidth and (TiT)^^ below Tc indicate that sample SI 
is also superconducting. 

First, we discuss the different behavior of the Knight 
shift K. in the superconducting state. When a supercon- 
ducting condensate consists of singlet Cooper pairs, a mag- 
netic field cannot polarize these paired electrons unless 
pairs are broken up. Consequently the spin susceptibility 
Xspin and thus K. vanishes [39] below Tc when T — 0. This 
may explain the behavior that we observe in sample S2, 
and that is reported for powder samples of LiFeAs [TOlf^ 
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so far. In contrast, the constant Knight shift across Tc in 
sample SI for H \\ ab exhibits a behavior that would be 
expected for a spin-triplet superconductor where pairing 
does not interfere with the magnetic response of the elec- 
trons and Xspin remains constant across Tc down to zero 
temperature. For small angles off H \\ ab, however, the 
Knight shift decreases below Tc, suggesting that a small 
component of the magnetic field along the c direction leads 
to a transition into a state with a superconducting or- 
der parameter of different symmetry in the orbital and/or 
spin sector. Note that the so far reported Knight shift 
measurements of LiFeAs [TOll^ have been performed in 
polycrystalline samples in which the broadening due to 
large second order quadrupole shift may hinder determin- 
ing the intrinsic spin shift if there is a strong angle depen- 
dence of the Knight shift. Another explanation would be 
that the nature of superconductivity may strongly depend 
on the effective doping level of the samples, even though 
LiFeAs is considered to be a stoichiometric compound. 
The differences in the doping level could be visualized by 
the slightly different quadrupole frequencies: those sam- 
ples with a high quadrupole frequency reveal the constant 
/C across Tc, whereas those samples with a lower vq, in- 
cluding the powder sample reported in ref. [20], show a 
more normal behavior. 

A measurement of the intrinsic spin susceptibility in 
the superconducting state via the Knight shift is often 
not straightforward and other reasons may lead to a con- 
stant K."'^. First of all, the correction for the large second 
order quadrupole shift could mask a decreasing Knight 
shift. However, we observe that the quadrupole frequency 
exhibits only an insignificant temperature dependence at 
low T, and the constant IC""^ is already recognizable in 
the raw data without quadrupole correction in Fig. 4 (a). 
Also demagnetization effects are negligible for || ab due 
to the plate-like shape of the crystals. Furthermore, de- 
magnetization effects reduce the actual field in the super- 
conductor and therefore should lead to a further decrease 
of the total shift K, below Tc- A simple heating effect of 
the sample by the radio frequency pulses can be excluded 
since we observed the decreasing Knight shift for small 
angles off H \\ ab and for sample S2 under similar condi- 
tions. Another reason for a constant /C in a singlet super- 
conductor could be an enhanced magnetic susceptibility 
via strong orbital magnetism |40j or via spin-orbit scat- 
tering in the presence of disorder |41| . Such a scenario 
appears rather unlikely. First, the strong angular depen- 
dence of the Knight shift in the superconducting state (see 
Fig. 5) is incompatible with a large non-spin contribution 
to /C"'', such as strong orbital magnetism. For the same 
reason the unchanged Knight shift through Tc cannot be 
explained by an impurity-based scenario. Moreover, our 
samples appear to be clean single crystals, even though 
tiny differences between the crystals have to exist. 

In LiFeAs the (7r,7r)-nesting and static antiferromag- 
netism are absent [S] . Investigating theoretically the mag- 
netic and pairing instabilities in an electronic model that 
incorporates the poor nesting properties and unusually 
shallow hole pockets of LiFeAs, Brydon and coworkers 



[15| find ferromagnetic fiuctuations to drive an instability 
toward spin-triplet p-wave superconductivity. These mag- 
netic fluctuations are related to LiFeAs being in the vicin- 
ity of nearly ferromagnetic, incommensurate, long-range 
spin ordered phases, the stability of which is governed by 
the detailed electronic density and interaction parameters. 
Such a spin-triplet scenario may be supported by the con- 
stant spin susceptibility across Tc observed for H \\ ab 
and the increasing linewidth with decreasing temperature 
in sample SI, although samples S2 and S3 clearly show 
a spin-singlet behavior. As we have discussed above, the 
experimental results in single crystals S1-S3 should result 
from intrinsic properties rather than extrinsic ones. This 
suggests that different phases are competing in LiFeAs 
where relatively small changes due to external fields, struc- 
ture or stoichiometry can cause transitions between un- 
conventional and more conventional superconducting ground 
states, and may also explain other contradicting experi- 
mental results in literature (see also Section [I}. 

Another anomalous behavior is the strange upturn of 
(TiT)~^ in sample SI and in the polycrystalline sample. 
Nakai et al. [5S] suggest four different possibilities for a 
similar upturn in Lao.srCao.iaFePO: (i) impurity contri- 
butions which can be excluded for both compounds, as ar- 
gued above. Furthermore, the relaxation curves are single 
exponential in LiFeAs, too, whereas impurities would lead 
to a multiexponential behavior of the relaxation curves, 
(ii) Vortex contributions to (TiT)^^ which can be ex- 
cluded for LiFeAs since the upturn occurs in zero mag- 
netic field (NQR). (iii) Slowing of magnetic fluctuations 
due to the opening of the superconducting gap. In this 
case, scattering of local magnetic moments with conduc- 
tion electrons in the normal state enhances the energy 
scale of the fluctuations. In the superconducting state the 
scattering is suppressed due to the formation of Cooper 
pairs, and the local magnetic moments may slow down. As 
argued by Nakai et al., such a behavior has not been re- 
ported yet. (iv) Collective modes of the spin-triplet pairs 
could give rise to novel spin dynamics in the SC state [38l 
142] . In particular, we note that the possibility (iv) may ex- 
plain the observed constant Knight shift in the very same 
sample SI. 

In order to reconcile the two different results even in 
the single crystals grown under the same conditions, we 
conjecture that the nature of spin fluctuations, which is 
generally thought to mediate Cooper pairs in non-BCS 
unconventional superconductors, could be either antifer- 
romagnetic or ferromagnetic depending on the proximity 
to a ferromagnetic instability in an extremely sensitive 
fashion. It is also interesting to note that there is a seem- 
ing trend that the EFG and thereby the NQR frequency 
i/Q are larger, when the Knight shift and (TiT)^^ exhibit 
the peculiar behaviors. This indicates a possible connec- 
tion between the doping level, which could be propor- 
tional to the EFG as in other iron pnictides, and a fer- 
romagnetic instability. Consistently, Co doping leads to 
a lower quadrupole frequency (EFG), and to a reduction 
of (TiT)^^ in the superconducting state, putting the Co 
doped sample further away from the ferromagnetic insta- 
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bility. This observation is also consistent with theory [111 
[15], where the sniaU hole pocket at the F point drives the 
ferromagnetic fluctuations. Regarding Co doping as elec- 
tron doping, the hole pocket should shrink upon doping, 
thereby weakening the ferromagnetic fluctuations |43| . 

5 Conclusion 

We have used NMR and NQR to investigate different 
LiFeAs samples. We find that tiny differences in the sto- 
ichiometry of the samples exist which lead to different 
normal state properties as well as to different supercon- 
ducting properties. A possibility to distinguish the sam- 
ples is the quadrupolc frequency, vq. Sample SI yields a 
slightly higher i/q than sample S2, whereas i'q of a Co 
doped sample is even below that of S2. At the same time, 
the Knight shift of sample SI is constant across Tc for 
H II ab. Although this may suggest spin triplet supercon- 
ductivity, we find that the Knight shift drops below Tc for 
small angles off the ab plane. On the other hand the Knight 
shift of sample S2 decreases in the SC state regardless of 
the angle which is compatible with standard singlet pair- 
ing of the Cooper pairs. The temperature dependence of 
the lincwidth varies also with the crystals. The linewidth 
of sample SI increases with decreasing temperature, indi- 
cating a growing of magnetic fluctuations, while that of 
sample S2 is temperature independent. Consistent with 
the observation of /C"'' = const, across Tc could be the 
increasing {TiT)~^ in the SC state of sample SI as well 
as of the polycrystalline sample which also yields a higher 
vq. The increase of {TiT)^^ below Tc could originate from 
non- vanishing spin degree of freedom in the SC state which 
give rise to novel spin dynamics. 
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